Correlations in the Sine-Gordon Model with Finite Soliton Density 
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We study the sine-Gordon (SG) model at finite densities of the topological charge and small 
SG interaction constant, related to the one- dimensional Hubbard model near half- filling. Using the 
modified WKB approach, we find that the spectrum of the Gaussian fluctuations around the classical 
solution reproduces the results of the Bethe ansatz studies. The modification of the collective coor- 
dinate method allows us to write down the action, free from infra-red divergencies. The behaviour 
of the density-type correlation functions is non-trivial and we demonstrate the existence of leading 
and sub-leading asymptotes. A consistent definition of the charge-raising operator is discussed. The 
superconducting-type correlations are shown to decrease slowly at small soliton densities, while the 
spectral weight of right (left) moving fermions is spread over neighboring "4/cf" harmonics. 
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I. INTRODUCTION 



In this paper we study the sine-Gordon (SG) model 
at finite densities of the topological charge and small SG 
interaction constant. The exactly solvable sine-Gordon 
model has a large number of applications in different 
subfields of condensed matter physics and statistical me- 
chanics. The particular case of the SG model of our 
interest below arose previosly in various studies of the 
one-dimensional Hubbard model close to half-filling, |l]] 
commensurate-incommensurate transition, g strongly 
correlated systems in higher dimensions , collective ex- 
citations in a multi-fluxon Josephson junction |4|, super- 
conducting S or quantum Hall B double-layer systems 
in parallel magnetic field and phase excitations in ferro- 
electric liquid crystals M. 

Theoretical investigations of the quantum SG model 
with the finite density of topological excitations (kinks) 
were undertaken by Haldane S, Caux and Tsvelik, || 
Papa and Tsvelik. pQ] These studies, employing the anal- 
ysis of the Bethe anzats equations, established the disper- 
sion of low-lying excitations and the so-called Luttinger 
parameter, defining the asymptotic behavior of the cor- 
relation functions. 

In the present work, we use an alternative method 
to investigate this problem, which was introduced by 
Dashen, Hasslacher and Neveu, [[llj and is usually called 
the modified WKB or semiclasscal approach. In this 
method, one first finds the classical solution to the equa- 
tion of motion and then analyzes the quantum fluctu- 
ations around this solution. Q It is known that this 
method gives excellent agreement with the available ex- 
act solutions of several models, and particularly the SG 
model. 

With the use of the modified WKB approach, we find 
that the spectrum of Gaussian fluctuations in our model 
gives the spectrum coinciding with the one obtained ear- 
lier by the Bethe anzats method. PjHl At the same time, 
applying the standard semiclassical treatment, we en- 



counter strong infra-red divergencies in higher orders of 
perturbation theory. The difficulties arise also with the 
definition of the so-called dual field associated with the 
charge-raising operator of the theory. It is known that a 
troublesome difficulty with the semiclassical approach is 
the existence of fluctuations bearing zero energy, which 
corresponds in our case to the translational symmetry 
of the problem. The treatment of these "zero modes" is 
best done by the collective coordinate method |l3| , whose 
straightforward application, though, does not cure our 
problem. It is shown below however that it is possible to 
use a somewhat modified collective coordinate method, 
to arrive at the action, free from the infra-red divergen- 
cies. Technically it is done by allowing the fluctuating 
field to enter the argument of the classical solution. The 
stable infra-red action is achieved at the expense of a 
more complicated form of the correlations. 

Proceeding this way, we obtain the leading and sub- 
leading asymptotic terms in the density-density corre- 
lation functions in agreement with earlier predictions by 
Haldane. J8J The definition of the charge-raising operator 
now becomes free from ambiguities and we can calculate 
the leading behavior of the superconducting-type corre- 
lations. Remarkably, the decay of the latter correlations 
is slow at small densities of kinks, an observation which 
confirms earlier results. 
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We demonstrate a non-zero quantum average of the 
four-fermion umklapp operator in the Hubbard model 
away from half-filling. Particularly, it results in the 
spectral density of the right-moving fermions being dis- 
tributed over the neighboring Akp harmonics of the den- 
sity. 

The remaining part of the paper is organized as fol- 
lows : We obtain and analyze the classical solution in 
Sec. |0|. The fluctuations around it are discussed in Sec. 
p| . The difficulties with the standard application of the 
semiclassical approach are outlined in Sec. IV. A way of 
circumventing these difficulties is discussed in Sec. M, we 
present the form of the correlation functions here. Con- 
cluding remarks are found in Sec. VI. The essence of the 



collective coordinate method is given in the Appendix. 



II. CLASSICAL SOLUTION 



We consider the Lagrangian density of the form 



£ = \&<, 



1 m 2 B 

^^-—(l-cos^ + H^d.. 



a 



2tt 



(1) 



on a line of length L. The SG interaction constant is 
assumed to be small, (3 <C 1, and m is the first breather 
mass. The term with the chemical potential H does 
not enter the equation of motion and defines the bound- 
ary conditions for the field </>. We write the topological 
charge density p — w^d x ((3(f)) and the average density 
p = Q J ' L. Using the well-known correspondence between 
the fermions and bosons in one spatial dimension, |1J] 
we associate p with the fermionic density. Particularly, 
cos/30 appears as an umklapp term in the charge sec- 
tor of the Hubbard model and p measures the deviation 
from half-filling, as 2np = Akp — 2ir. The total charge Q 
is given by 



pdx 



h fi +ti> 



(2) 



which defines the boundary conditions for (f>. 

In what follows, we employ the quasiclassical method 
of analysis of the quantum system defined by (Q). This 
method arises most naturally in path integral formulation 
of the problem and basically corresponds to the steepest 
descent method. At a first step, one finds the classical 
configuration of the field, delivering the extremum to the 
exponentiated action, and then analyzes the spectrum 
of fluctuations around it. The general criterion for the 
applicability of the method is the formally large value 
of the classical action; it was shown, however, that the 
quasiclassics gives exact results for a particular case of 
sine-Gordon model with a zero density of the topological 
charge. Jll[ | 

For later convenience we make a shift /?</> — > /3<f> + n 
writing the potential term in the form 



U($) = Jl(l + cos/ty) 



(3) 



As a first step we seek a static classical solution </>o by 
variating C in <5</>. We have 



d 2 x fo - U'((f> Q ) = 



(4) 



Multiplying it by d x 4>o and integrating over x we obtain 

(c^o) 2 = 2U(4> ) + C 

with some constant of integration C, which is integrated 
again to give 



dy 



^2U{y)+C 



Letting C 



Am' \-k 



'■£— with yet undetermined k, we get 



"P -- *? 



X = — F(8d>o/2, k) + const, 

771 



with incomplete elliptic integral F. [[L5| Inverting the last 
equality, we have 



0o = 2/3 1 am(a;777/fc, k) 



(5) 



where am(i, k) the Jacobi amplitude function fllq] with 
the elliptic index k. This index is determined by the 
above total variation of the field (/3</>/2) which results in 
the equation p6| 



£ = kK(k) 
2p 



(6) 



with the complete elliptic integral K . Henceforth we will 
omit the index fcasa second argument of the elliptic func- 
tions and use the conventions k\ = \fl — k 2 , K(k) = K, 
K{k x ) = K', E{k) = E. From (g) we have 



fci ~ 4e 



-rn/2p 



P <C 777, 



1, p > 777. 



(7) 
(8) 



Introducing the soliton mass M s = 8m/ p 2 and rescal- 
ing the field H = M s h we write the energy of this classic 
solution, £q = —C[4>o] as follows 



P 2 £ _2E- k\K 2h 
~ kK 



2m 2 L 



k 2 K 



(9) 



The index k, defining the density (number of kinks), 
should be chosen in order to minimize the energy £q. 
Differentiating m) by fc, we find that the energy attains 



its minimum at 



E/k = h 



and this minimum has a form 



So 



2m 2 L k\ 
~~f3 2 ~~k 2 



(10) 



(11) 



The last expression coincides, up to the sign, with the 
pressure p = —d£o/dL\Q. The enthalpy of the system is 
given by ^^. 



A. Number of Kinks 



III. FLUCTUATIONS AROUND THE CLASSICAL 
SOLUTION 



At smaller fields, h < 1, there are no solutions to Eq. 
(|l0|). The critical value of the field for the appearance 
of the "kink condensate" , H cr = M s was found earlier in 
|^0[ . At fields, slightly exceeding this critical value, we 
have for the density (m/p)e~ m ^ ~ (h— l)/4. Upon fur- 
ther increase of the field, we have approximately linear 
relation between the field and the density of kinks : 



in I In 

j2 



AH r 



H -H c 

.2\ 



H ~ H r 



Hp'/lAar*), H > H e 



(12) 
(13) 



Numerically, the asymptotic expression for large kink 
densities (13) holds with a good accuracy already at 

Pt > 2H rr . 



B. Susceptibility 



Consider now the susceptibility defined as \ = 
—L~ 1 d 2 £o/dH 2 . A simple calculation gives 



f3 2 E 
16 k 2 K 3 

P 2 
8(h- 1) V" h- 1 
(/3/2tt) 2 , h > 1 



In - 



(14) 

(15) 
(16) 



Our expressions for the kink density p and susceptibil- 
ity x are identical with those obtained earlier [|l0| in the 
large-field limit, eqs.(l^), (|l6|). At the fields H ~ H cr our 
expressions for p and \ coincide with their counterparts 
in [1(|] up to an overall factor 7r/2. 

Near its minimum at p = p, the energy can be ex- 
panded as follows 



£o 



2.Y 



(P-Pf 



(17) 



or, 1/x — ftVN, in the notations by Haldane. |ra| 

Another important parameter vj, according to Hal- 
dane, is given by the coefficient in the total Hamiltonian 
H = vjP 2 /(2irLp 2 ), at small values of the total field mo- 
mentum P = J Q dx<p(p'. It can be easily shown, that 
H = P 2 /2A4 t ot with the mass of the kinks' condensate 
Mtot = Jodx(d x (f>of = LE/K{2m/(3k) 2 = QH. As 
a result, one finds vj — np/H — nf3 2 /(16KE). There- 
fore, one expects || the Fermi velocity vf — y/vNVj = 
kiK/E and the Luttinger parameter /C = s/vj/vn = 
n(3 2 /(16kiK 2 ), both expectations verified below. 



A. Lame Equation 

Now we consider the fluctuations around the classical 
solution fa. We write <fi = cf>o + r\ and expand the La- 
grangian into Taylor series pjfl 



C = C[<p ] + d 

£i = \{dtV? \(d x v) 2 + l -m 2 {cos[3M v 2 
+higher orders in r\ 



(18) 
(19) 



The term linear-in-r; drops out from this expression and 
the higher orders in r\ contain additional powers of /3, 
and we will neglect them for the moment. We look for 
the solutions r\ in the form r](x,t) = J2 c i(^) ? 7i( a; ); obey- 
ing cyclic boundary conditions. Integrating by parts, we 
write 



A = i,(- 



-d 2 + d 2 + m 2 cos /3(j> ) V 



Next, we seek the "normal modes", satisfying 
(-d 2 - m 2 cos @ '(j>o)r]i(x) = uJ 2 rji{x) 



(20) 



(21) 



The dynamics of these normal modes is simple, c(t) 



?±iujt. 



Introducing the variable z = xm/k we rewrite 



Eq. (H) in the form §] 



dz- 



: r) = (A + 2k 2 sn 2 z)r] 



(22) 



with A = —k 2 (u> 2 /m 2 + 1) and sn(x) Jacobian elliptic 
function. This is the Jacobi form of the Lame equa- 
tion and the solutions to it are described in the liter- 
ature . (Ref. p7| , chapter 23.71) These solutions are 
parametrized by an index a and are explicitly written as 



. . H(z — a) 
V{z,a) = exp(zZ(a)) 



(23) 



with Jacobi functions H(u) = l di(iru/(2K)), O(u) = 
r d4(TTu/(2K)) and Z(u) = dln<d(u)/du; the nome of the 
i}— functions q = exp(— nK'/K). jl5| The structure of 
( p3| ) shows that H/8 is the modulating Bloch function, 
and Z corresponds to a wave vector. The energy u> is 
given by 



uj = ±(m/k) dn(a) 



(24) 



The second independent solution for n with the same en- 
ergy ( p4] ) is obtained from ( p3| ) by changing a — *• — a 
which corresponds to the symmetry of eqs. (p2|), (|2J 
with respect to reflection z — > — z. 



B. Properties of 77(2) 

One can show that the whole family of linearly in- 
dependent solutions (p3|), which possess the real- valued 
energies and are periodic on a length of the chain, is 
exhausted by the values of a belonging to segments 
{K - 2iK',K] and \-2iK',0) in the complex plane. We 
will refer to these segments as to segment I and II, respec- 
tively. The energies ( J24[) and the normalized solutions 
(B3T) are doubly periodic in a with the periods 2K, 2iK'. 

The solutions ( |23| ) are quasi-periodic in 2K, and 



rj(z + 2K, a) = -r)(z, a) exp(2KZ(a)) 



(25) 



Introducing the Floquet index v by e 2lvK = —e 2KZ ( a \ 
we have 



7r mi 



(26) 



with integer n. Returning to the original variable x 
zk/m, we find the allowable wave vectors in the form 



q = 2irp(i(K/ir)Z(a) + 1/2) = 2nn/L 



(27) 



Consider first the variation of the index a on the segment 
I, the corresponding n in (E7J) in the range 



- Q/2 < n < Q/2, 



(28) 



which means Q allowed values of n. The energies ( p4| ) of 
solutions r)(z, a), with a from the segment I, lie between 



and 



lo = 0, q = at a = K± iK' , 



lo = u)\ — mk\/k, q = np at a = K. 



(29) 



This part of the spectrum is interpreted as a band, orig- 
inated due to hybridization of the bound states related 
to individual kinks. 

In order to see this, we note that for one kink, Q = 1, 
the only value n — is allowed in the segment I. It is 
the quantized bound state with to = for the single-kink 
solution of the sine-Gordon equation. In the limit of low 
density p <C to, the energies of the lower band are expo- 
nentially small, 



4me-' m ^P sin' 5 



2p> 



(30) 



which is the dispersion for a system of weakly coupled 
harmonic oscillators. 

Another branch of the spectrum, parametrized by a 
from the segment II, has one singular point at a = —iK', 
where both Z(a) and dn(a) have simple poles. This point 
corresponds to q — > ±00 and uj — > 00. 

The energy of this band has the minima 

Lu = uj2 = rn/k at \q\ = ttj5, a = —2iK',0. (31) 



The dispersion at p/m = 0.26 is shown in Fig. [j]a in 
the extended Brillouin zone scheme. The same dispersion 
shown in the reduced Brillouin zone scheme is depicted 
in Fig. 0b. 




-0.5 0.5 

q/7lp 

FIG. 1. The dispersion in the (a) extended and (b) reduced 
Brillouin zone scheme shown at k = 0.87 and p/m = 0.26. 

We see that there is an energy gap LO2 — lo\ between 
the states in the lower and upper band. It is interesting 
to note a following relation 

22 2 

LU 2 — f^l = TO . 

In the low-density limit, p ^C to, we write a = ia' ', then 
q — > —to tan a', u> — > to/ cos a' so that ui — y^m 2 + q 2 . 

At this point it is worthwhile to calculate the Fermi 
velocity vf for low- lying excitations in the large-Q limit. 
Expanding u and q near the point ao = K — iK' of sector 
I, we have 



. fel, s 

lo ~ —tm—(a — ao), 
k 



-i2pE(a - a ), (32) 



and 



vp = lo J q = k\K / E, 



'(/i-l)ln 



h-1 



1, 



l -\M) ' ft>>1 ' 



(33) 

(34) 

(35) 



in agreement with the above estimate and the results in 
|L0|]. As a result, the Luttinger parameter, K, — irx v F, is 
found in the form : 



which means 



K 



7T/3 2 



1 



16 kiK 2 ' 



(36) 



which shows that /C decreases with increasing the density, 
and attains its limiting value /3 2 /(4_7r) at p — * oo. 

The form of the eigenfunctions (|2^) is somewhat sim- 
plified at the special points a = K — iK 1 , K, 0, where one 
has (cf. ©, fl|)) 



rjo(x) = (LE/Ky 1/2 dn(xm/k), 
r)i(x) oc cn(xm/k), a = K 
772(2;) oc sn(xm/k), a = 



K-iK' 



(37) 
(38) 

(39) 



with the normalization factor in 770(2;) made explicit. 
Henceforth, we will call the eigenfunction (g7|) the zero 
mode. 

Analyzing the lowest-energy states with a ~ K — iK' , 
one can find a following expression, valid in the leading 
order in q/np : 



rj q (x) = rj (x) exp -iqx - i 



qn 



■d 



AKEp d 3 (nxp)J 



(40) 



The second term in the exponential is the phase of the 
Bloch function. It is 1/p— periodic in x, vanishing at 
x = 0,±p~ 1 ,±2p- 1 ,.... 

For completeness, wc provide here an approximate ex- 
pression for the eigenfunctions of the upper band in the 
low-density limit, p <C m and k — ► 1. Near the position 
of nth kink, it is given by 

vi x )u~m oc (tanh 777,(2; — pn) — itana) 

x exp(z(x — pn)m tana — inq/p), (41) 

q = — mtana + lap + np, (42) 

in accordance with flial . 



Concluding this subsection, we note that eqs. (33) 
(pq), (|40|) , ( |4l| ) extend the analysis of fluctuations per- 
formed by previous authors 



C. Quantization of Fluctuations 

The differential operator in (E2) is self-adjoint and, 
consequently, the eigenfunctions form a orthonormal set, 



dxr) n (x)rj m (x) = S n 



(43) 



There is a well-known property that a continuous func- 
tion F(x) (with the same boundary conditions) can be 
expanded as a generalized Fourier series 



F ( x )=y2vn(x) dyr) n (y)F(y) 
Jo 



(44) 



^2vn(x)Vn(y) =6{x-y) 



(45) 



Now we may introduce the quantized field, defining 



${x,t) = Y, -4= (Vn(x)e^ t bl + r&(x)e-*>» t b n ) 



(46) 



The Bose operators b n satisfy the commutation relations 



[bniblj =5 n;T 



and therefore for equal-time fields 



(47) 



Mx),dt$(y)] = lJ2^n(xWn(y)+V:( x )Vn(y)) (48) 



iS(x - y) 



(49) 



Strictly speaking, the zero mode (B7p with lo = should 
be excluded from the sum (M), which leads to the ex- 



tra term —ii]o(x)rio(y) ~ Q~ l — > in (|48|). The reason 
for it is non-independent character of the mode 770 in the 
decription of the field configuration, 4>o(x), {rin(x)}^Lo- 
Indeed, 770, bearing the zero energy, corresponds to the 
translation of the classical solution <fio(x) — ► (f>(x + X), 
while 770(2;) ~ d x 4>o(x). As a result, the consistent treat- 
ment will consist of the functions <po{x + X),{7] n (x + 
X)}^ =1 , with 770 replaced by the displacement X of the 
system as whole. The quantity X should be regarded 
as a dynamical variable in the correct description of the 
quantum Hamiltonian. |13| 

We outline this collective coordinate method in Ap- 
pendix and show, that the corrections to the above sim- 
pler formula are negligible in the limit Q — ► 00. The 
only thing, which we should borrow from the collective 
coordinate method, is the necessity to integrate over the 
vacuum location X instead of inclusion of 770 . 

It seems that one now has everything to compute vari- 
ous correlation functions, describing the quantum fluctu- 
ations of the system. Given a set consisting of the clas- 
sical solution <pQ and the normal modes around it, one 
calculates the quantum averages with the simple rules of 
bosonic algebra and then integrates over the initial po- 
sition X of the soliton lattice. This program however 
reveals many difficulties described in the next section. 



IV. DIFFICULTIES IN THE APPROACH 

A. Density- type correlations 

Consider the fermionic density p = d x (3(j)/2ii. The clas- 
sical contribution po = d x l34>o/2ir = pw(x) where 



The lattice of kinks gives rise to the notion of the Bril- 
louin zone and to the infinite set {?y„ '} of the eigen- 
states with non-zero energy and zero wave vector in 
the reduced Brillouin zone scheme; t]q = r]o. A con- 



stant is not an eigenstate of (21), and therefore rj n 



L 1 J Q dxrjn (x) 



, . IK , m(x + X) 

w(x) = (in ; 

7T K 



(50) 



The averaging over the kinks position gives a simple re- 
sult here, 

{p{x))x = p(w(x)) x = P- 

At the same time, the classical contribution to 
the pairwise density correlator is (p{x)p(y))x = 
p 2 (w(x)w(y)} x 7^ p 2 ■ These non-trivial correlations of 
density at largest distances correspond to Wigner crys- 
tallization of the charge Q. This long-range ordering of 
the charge in one spatial dimension is evidently erroneous 
and should eventually be removed from the theory. 

The quantum corrections to the above correlator in the 
long-distance limit read as 



r,- 1 / 2 . Further, if the wave vector 
q n of the eigenstate r\ n is close to zero, a correspond- 
ing Fourier transform rj n {q) has a pole fj n /(q — q n ). The 
quantity fj n decreases when increasing the energy and 
this decrease is slower at larger k. In particular, one can 



show for the above set {rj 

while ?7o 



(0) 



} that XT=0 (Vn)) = L- 



tt/(2VLKE). 

The vacuum expectation value of the product of 9 fields 
is then 



(0(x + y)9(y)) 



E V F _ 2 
— ?7 cosg:r 

q>0 q 



Ecosq x 
q' 



2q 



12 



£< 



'n.q'Vn' 



(53) 



(p(x)p(y)) q = 



K d 2 

2ir 2 dxdy 



where we have used @, @, @, (j4fi 

The CDW-type correlator is similarly estimated as 



As usual, the constant (divergent) term should be added 
here in order to analyze the dependence of the correla- 
tions on x. The first term in (p3), resulted from the 
lowest energy modes, contributes to the logarithm. In 
the second term, with both signs of q' available, we note 
(w(x)w(y))x In \x — y\ (51) that linear-in-g' corrections to u) n ,q' vanish and the sum 

over q' gives a contribution linear in |ar|. The sum over n 
in (p3) converges rapidly when elliptic index k is not too 
close to unity, i.e. at p ^ m. For an estimate, it is suffi- 



(e 



0<t>( x ) e -W(j/)\ 



(exp(-2/Cln(x - y)w(x)w(y))) x 

(52) 



Evidently, the modulation w(x)w(y) in (|52| ), being expo- 
nentiated is more pronounced than in Q51|). 



B. Dual Field 

The so-called dual field 6(x,t) used in the construction 
of the fermion operator, is usually defined as 



d x 9 = d t <t>, d t 9 = d x 4>. 



>iVl 



p(L~ 



V 2 o) 



cient to write J2 n >i w ",g"7n ~ P2~2 n _ 
As a result, we have 

-p\x\(l-TT 2 /(4:KE))0{l) (54) 

The linear in |x| term is present unless k — > (raj p — > 0), 
when one recovers the usual expression, (9(x + y)9(y)) ~ 
— ln|x|. As a result, the superconducting correlations 
^ e 47ri/3 e e -Airip 9^ decay exponentially at large dis- 
tances. 

We see that the above definitions of the dual field are 
inconsistent in the basis r\ n and lead to an exponential 
decrease of the SC correlations in the gapless situation. 



These relations imply d 2 <\> — d 2 <fi, which is not consistent 
with the linearized theory equation ( |2l|) (the discussion 
of the zero mode contribution to 9 is postponed until 
the next section) . Particularly, the usual construction of 

the charge-raising operator O ~ exp 2iri(3~ 1 J Xl (dtd x + 
dxdt)4> is inappropriate in the basis {ry n }, being depen- 
dent on the contour of integration in O. 

One may ignore for a moment this discrepancy, using 
only a first equation d x 9 = dt<j>. This is sufficient for 
the calculation of the instantaneous correlation function 
(exp(47ri/3^ 1 6'(a;))exp(-47ri/3~ 1 6'(y))). Making this cal- 
culation, we meet another troublesome property. 



C. Perturbation Theory 

Consider now the role of the terms of the interaction 
in (119) , which have been discarded. These terms, con- 
taining additional /3s, produce divergences in the per- 
turbation theory. Resumming the tadpole sequence of 
diagrams, Fig. 0a, one obtains the mass renormalization 
of the form 

if \x)m 2 cos I3(j) {x) -» 7 ? 2 (a;)TO 2 cos0o(x)e-' c,i,(a:)2lnA 

with A ~ Q. This renormalization depends on x and, 
therefore, changes the shape of the potential cos/30o- 



This essential modification is accompanied by stronger 
than logarithmic infra-red divergencies in diagrams, not 
contained in the tadpole series, e.g., that shown in Fig. 



|b 



H 




(a) 

+ . . 



FIG. 2. (a) A tadpole sequence of diagrams, contributing 
to position-dependent mass renormalization. (b) A stronger 
divergent diagram not contained in this sequence. 



V. A PROPOSED SOLUTION 

A. Lagrangian 

A key to resolve all the above discrepancies is found as 
follows. We wrote the fluctuations around the classical 
solution as = 0o + V and performed the Taylor expan- 
sion (|19|). Instead of that, we can write the fluctuations 
in a form 



4> = 2pr\ 



2K[xp + ^ 



(55) 



The formula (|35|) has a remarkable property, that shift- 
ing fj — > fj + 2ir/f3, we increase accordingly <fi — > (j) + 2n/f3. 
This may serve as a justification of the scaling factor be- 
fore fj in (p5[). At p > m one has <f> — > 2irp(3~ l x + fj, 
hence (q5|) and the previous writing cf> = tfio + r\ coincide 
in this limiting case. Some calculation shows, that the 
perturbation theory in fj can be obtained from the exact 
form of the full Lagrangian 



C 



j3v_ 
2np 



(dtv) 2 - (d x fj) 2 



(56) 



The field theory decribed by ( pq ) is free from infra-red 
divergencies for two reasons. The first is the presence 
of the field derivatives in each term of the perturbation 
series. Another reason is a particular form of the form- 
factors of the interaction, when nth vertex is obtained by 
(n — 2) differentiations of the restricted function w 2 (x). 
The Gaussian action is given by 



(57) 



(58) 



C = \w 2 (x) ((d t f,) 2 - (d x fj) 2 ) . 
with the corresponding equation of motion 
d 2 fj = w~ 2 d x w 2 d x fj 



Our aim is now to show that the spectrum of (pq) co- 
incides with one of the initial problem ( |2l| ) , and to find 
the connection between the fields fj and r\. 

The equivalence of the spectra stems from the possi- 
bility to rewrite the Lame equation in the form : 



u! 2 rj(x) = —w l (x)d x w 2 {x)d x w (x)r)(x) 



(59) 



which implies that the Lagrangian is 



1 



C = g {{dtv? ~ (wd.w-^) 2 ) 



(60) 



Evidently, eq. (|60|) coincides with (pTf ) after the substi- 
tution 



fj(x,t) — w l (x)rj{x,t). 



(61) 



As a result, the eigenfunctions to (p7|) are normalized 
with a weight 



dxw 2 {x)ij a {x)f]f3(x) = 5 a/ 3, 



(62) 



and ( pl| ) leads to the following form of the correlations 
at large distances 

1 (fj(x, t)fj(y, 0)) ~ -K In | (x - y) 2 - v 2 F t 2 \ (63) 

Therefore, we find the description of the fluctuations, 
possessing the same spectrum in the Gaussian approx- 
imation, but free from divergencies in higher orders of 
interaction. This choice of appropriate variables should 
be hence considered as an effective partial resummation 
of perturbation series. 



B. Density-density correlations 

The adoption of the definition ( p5| ) makes the calcu- 
lation of the correlation functions less trivial. Consider 
first the density-density correlations (p(x)p(y)). Using 
the series representation Ea] for the function dn(x) we 
write 



(S+^&S 



E ' 



in(2-irpx-\-f3r)) 

coshnr 



(64) 



with r = ttK' /K. Note that the previous kinks' lat- 
tice displacement X corresponds to the new zero mode 
rjo = const, eq.(|6l|). The integration over fjo gives evi- 
dently (p) = p, while it selects the terms with equal ns 
in the pairwise density correlator, 



(p(x)p(y)) q = p 2 



K 



J2 Cn 



2ir 2 (x-y) 2 
cos2mrp(x — y) 



\x-y\ 



2n 2 K 



(65) 
(66) 



with omitted faster decaying terms. We see that the 
weight factors w(x) in ( |5l| ) are transformed into the de- 
caying Akp = 2ixp harmonics of the density, in line with 
Ref. i|. 



C. CDW-type correlations and umklapp ordering 

Keeping in mind the Hubbard model near half-filling, 
we associate the exponent e^^l 2 with the "2fci?" density 
wave operator, whereas e*^ stands for the "4/cf" umk- 
lapp process. 

The "2fcF" correlations are estimated as follows. We 
use (pq) and the series representation 



\{2Kz) 



r (l) 



V^ q(1) e 2wiz(n+l/2) 

71— — OG 

ne r(n+l/2) 

kKsmh(2n + l)r' 



(67) 



to obtain 



( e i#<«z) e -if<«0)\ 



E [<W 



n— — OQ 



p 2iir px (n+1 / 2) 
1^12^(71+1/2)2 



(68) 



(e 



i|0(x)_i|0(O)\ 



00 _2wrpx(n+l/2) 

V^ r (i) r (i) e rfiq \ 



n— — oo 



• | x |2K;(b+1/2) 2 



The prefactors of the exponents here are not universal, 
depending on the cutoff. What matters here, is the ratio 
of the amplitudes for the forward-going and backward- 
going waves with the same absolute value of (n + 1/2). 
The leading asymptotes are defined by n = 0, — 1 ; noting 
that C^l = -qC^\ with q = e-* K '/ K , we write 



Jf3<t>(x)/2 e -ip(f>(0)/2\ ^ r e iirpx + ~2 e ~iTr P x\ 11 -/C/2_ fjQ\ 

( e iM*)/* e *M°)/2) „ -2qcos(npx)\x\-^ 2 , (71) 



Particularly the above formulas indicate the following. 
The lowest harmonic of the CDW order parameter field, 
ex sin(/30/2) = sn[2K(px + f3fj/2n)] has an amplitude 
n/(kK sinhr/2) which remains finite in the total range 
of variation og p. The similar quantity for the SDW 
field, oc cos(f34>/2) = cn[2K(px + (3f]/2n)] 1 has an am- 
plitude it / (kK coshr/2) which vanishes as p/m at the 
half-filling, p = 0. 

Using ( |67| ) we write the umklapp order parameter field 
in the form 



cos 



(2) in(2irpx+Pfj) 



= E <$ 

n=—oc 

2C^ = E^„ 1 Cl 1 1 ) + (n-,-n). 



(72) 
(73) 



We observe the appearance of the constant term, 

(2) 

Cq , corresponding to the fact that the average value 
^iPMx)^ — {cos f3(j)(x)) — Cq ' is not zero at finite den- 
sity of kinks. Indeed one has 



C™ = {-2E + K{l + k{))/{k 2 K), 
~ 1 — Ap/m, p/m < 0.2 
~ m 2 /(8n 2 p 2 ), p/m > 0.2 



(74) 
(75) 
(76) 



(2) 

It should be noted that the non-zero value of Cq arises 
without "locking" of the field <f>. We remind that in the 
absence of kinks the field <j> is "locked" to one of the 
discrete values </>„ = 2nn/P, and it results in the non- 
zero value of the "2fcp" density wave (e I/3< ^ 2 ) = ±1. In 
the presence of kinks <f> connects different vacua <p n , and 
{ e l P<fr/ 2 ) — 0, while (e 4 ' 3 *) ^ is decreasing function of 
the kinks density. 

The quantum average of two umklapp exponents at 
long distances has a form 



(cos /3<j>(x) cos /3</>(0)) - [C, 



t cos 2irpxn) 



-(2)l2 , v^ ' ^P x '< 

i + Z^ C " \r\2Kn 2 
n=l 



(77) 



(2) 

We see that in the absence of Cq the dimensionality 
of the operator cos/3</) would be JC. As K, increases 
with decrease of p, one might think that umklapp op- 
erator cos (3(f) becomes less relevant, pi However, this 

increase of K. is accompanied by the increase of Cq , and 
(cos (3(f)) ~ 1 at p ~ m. 



D. Dual field 

Let us now discuss the modification of the definition 
of the dual field 9(x, t). The basic property of the charge 
raising operator O ~ exp27ri/3 -1 # in the path integral 
representation is that the contour of integration ( "Dirac 
string" ) in the definition of 9 introduces a discontinuity 
in the field </>, so that the values of the field (3(j> differ by 
2ir across this contour. ]19[ In view of the above property 
(3(f) — ► (3(f) + 2-7T as (3fj — ► (3fj + 2ir, we see that it suffices 
to require the charge-raising property for the field fj. 

The canonical momentum for the field f] is given by 
7r = dC/d{dtf}) — w 2 dtfj. The quantization condition 
reads as [rj(x),ir(y)] — iS(x — y). 

Consider the definition 



dx' tt{x') 



(78) 



for equal-time fields. Differentiating the last equality by 
x and then by t we have from (J^) dtd x 9 — d x w 2 d x rj, 
hence the definitions 



d x 9 = w 2 d t fj, d t 9 = w 2 d x f], 



(79) 



are consistent with the equation of motion. The charge- 
raising operator, O ~ exp27ri/3~ 1 J w 2 (dtd x +dxd t )fj, 
is now uniquely defined, independent of the contour of in- 
tegration. Further, one has [fj(x),6(y)] = M?(y — a;) and 
therefore we have ex.p(2nif3~ 1 9(x)) f (Tj(y)) — f(f)(y) + 
2irf3- 1 $(x-y)) exp(2nif3' 1 9(x)) for any f(z). This prop- 
erty allows one to define the fermion operator, see below. 
A counterpart 9q of ryo = const cannot be obtained 
from ([79), which shows again a particular role of the 



zero mode. One of the possible ways to introduce the 
quantized quantity 8q, relevant to our discussion in Ap- 
pendix, is found in p0[ . One can write these zero-energy 
components as follows : 



m = 



Jt+(q-Q)Wi(x) 



8LEK 

TT 2 f3 

8LEK 



.2d_ 
2n X L + l p~~dJ 



(80) 



JW 2 {x) + (q - Q)t 



P 



2n X L 



2vr dq y ' 



with J = P/itp the persistent current attaining inte- 
ger values, the functions W\(x) — f x w~ 2 (y)dy and 
Wb(aO = f X w 2 {y)dy are expressed through the elliptic 
functions Nd(x) and Dn(x), respectively. A canonical 
conjugate to J is irpX, so the last term in (|80| ) is sim- 
ply the previous shift in the initial kink position. The 



charge variable q was assumed to be Q in Sections III 



[V and we let it fluctuate now around its equilibrium 
value. A choice of the factors before J and q variables is 
dictated by the condition J Q dxw 2 [{d t f]o) 2 + (d x f)o) 2 = 
{q-Q) 2 X~ 1 L~ 1 +P 2 Mz},ci. (||1). A particular form of 
the function W\(x) in (pfl) is not coincidental here. Sat- 
isfying the equation w~^d x w 2 d x Wi(x) = 0, the function 
W\ is another (increasing) zero mode for eq. (pq), which 
corresponds to the increase of Q in the classical solution 
w{x)W\{x) tx d(j>o/dk ; we have (3fjo | x=0 = 2n(q — Q). 

Let us discuss now the correlations of the dual field 8. 
To do it, consider first its equation of motion, as follows 
from © : 



d 2 8 = w 2 d x w~ 2 dj 



(82) 



Recalling the property of the Jacobi function dn(s) , we 
write 



w[x - 



"72) 



C\W 



- l (x), 



C\ 



AK 2 k x 



(83) 



Therefore the differential operator in the r.h.s. of (|82|) co- 
incides with one in (pq), when shifted by its half-period, 
x -> x + p' 1 /2. 

For the infra-red Gaussian action, given some partic- 
ular solution fj a to (pq), we obtain its dual counterpart 
8 a with the same energy and Floquet index. In view 
of the completeness of the set fj a , we have a relation 
9a(x) — c a fj a (x + p _1 /2). The value of the constant c a 
is elucidated from the following sequence of equalities 



c a d t fj a (x + p 1 /2) 



d t a {x) = w (x)d x 7) a (x) 
e 2ivK w 2 (x)d x fj a (x + p- 1 ) (84) 
e 2wK w 2 (x)c~ 1 dJ a {x + p- 1 /2) 
e 2ivK c\Cd^ a {x- 



p-72) 



hence 



±cie 



ivK 



(85) 



One can show, that the sign here is plus for a € 
(K — iK',K) and a € (—iK',§) and minus for a S 



(K - 2iK',K - iK') and a e (-2iK',-iK'). In the 
large- p limit the eigenfunctions fj a (x) are usual plane 
waves, ~ e -iq<*x ^ anc j we re cover the usual expression, 
@a(%) = sign{q a )fj a {x) . The sign in ( p5[) is irrelevant for 
calculating the expectation value of the product of two 
dual fields and we have 



0(x)9(y)) = c 2 (fj(x)fj(y)). 



(86) 



Noting from (B6| 
ately obtain 



that 47rci = f3 2 /IC, we immedi- 



/ e 4 7 r 4 /3- 1 9( a; ) e -47ri/3- 1 e(y)\ ^ i _ y\'V^_ 



(87) 



m is ob- 

\2V ' 



The dynamics of the correlations (Bq) , (t 
tained by replacing |x — y| — > {{x — yf — v^Xti — tz) 2 ) 1 / 2 - 
This simple temporal dependence of the correlations is 
valid as long as we consider the long-time, long-distance 
behavior, x ~ vpt 3> p _1 . At shorter scales, the compli- 
cated structure of the wave functions and of the disper- 
sion comes into play and the Lorentz form of the corre- 
lations is lost. 



E. Fermionic correlations 



We now discuss the Hubbard model and regard cos 0(f) 
as four-fermion umklapp term. In this case, in addition 
to the "charge" field discussed insofar, one also consid- 
ers the spin degree of freedom described by the field 4> s 
and its dual one S . We write the right and left fields as 
the linear combinations 



Ra - 


■S*-7»-\ 


w 


>La = 




/I 



with the spin projection a = ±1. The operators 



* 



Ra 



Xcre 



i<t>Ra 



i>. 



La 



Xae 



-i<pLa 



(89) 



(90) 



may be identified with the right and left fermion, respec- 
tively. Let us discuss this point in more detail. 

One can easily check, that given the usual re- 
lation [cf) s (x), s (y)] = id(y — x), the operators 

^Ra-> V'-Ro-, V'lo-' ^ l<7 with the same a anticommute. The 
anticommutation of these operators with different projec- 
tion of spin is ensured by the Majorana fermionic vari- 
ables Xa- @1 These variables obey X+X+ — X-X- = L 
while x+X- = -X-X+- 

A subtler question concerns the relation of the con- 
struction (M) to the initial fermionic variables in the 
Hubbard model. In the bosonization procedure, the ini- 
tial fermions were written similarly to (|9C|), with the dual 
field 8 in ( p8[ ) instead of our 8. The oscillator represen- 
tation (Eq) for the initial right (left) bosons had simple 



form, which is achieved within our formalism for p ^ 
at m = 0. It is not at once clear that the complicated 
expressions ([35]), (179), ( |23| ) in the presence of umklapp 
term, m ^ are related to right and left moving parti- 
cles. We state here that the label "right" of the boson 
field (jjRc simply indicates the species of fermion, which 
is constructed this way. As we will see shortly, umk- 
lapp interaction produces a partial redistribution of the 
fermionic spectral weight over the 4fcF-components of the 
density, hence initially right-moving particle has now left- 
going components. 

In short, our fermions (^0|) anticommute, have the 
same initial field <j>, and explicitly acquire the conven- 
tional form at 77i = 0. Therefore they are the fermions, 
which one starts with during the bosonization of the Hub- 
bard model. 

The spin fields of the fermions are factorized, and for 
the charge part we have 



iPr,l(x) oc exp(±7-</> - i-rO) 



-27i-i/3 -1 9 



E <$\ 



±i(2Trpx+l3fj)(n+l/4) 






(91) 



(92) 



(93) 



cf . ( p7J ) , ( [73] ) . The fermionic correlations are found in the 
form 



(4a(M)lM0)) 



(x-v,t)Va 

OO 

E 



(94) 



c it e 2iripx(n+l/4) 



n— — oo 



In 



{x - v F t) 2n + 1 / 2 {x 2 - v 2 F t 2 )^ 
((2n + l/2)/C - l) 2 



4/C 



(95) 



In the case of free fermions (m = 0, 1 = 8ir) one has 
K, = 2, and the anomalous dimension of the fermion 
vanishes, 70 = 0. We observe that at special values 
of the Luttinger parameter, K./2 = 1/3, 1/5, 1/7, . . ., 
certain j n vanishes as well, and the fermionic correla- 
tions reveal purely chiral (sublcading) component. For 
K = 2/(277 + 1), it is of the form 

„±7Tipa;(n+l/2) 

(iTM) n+1/2 (i-w»t) 1/2 ^ 96 ^ 

with upper (lower) sign for even (odd) 77. The leading 
asymptote in (|94| ) is given by a term with n = and 
hence the scaling dimension of the fermion operator has 
its usual value dim{ipR a ] = \(^JJC/2 + ^2/JC) 2 . 



F. Low-density limit 

We have obtained the critical exponents for the cor- 
relation functions which depend on the Luttinger pa- 
rameter K.. Some inspection of the formula (pq) reveals 



that at small densities the Fermi velocity is exponen- 
tially small, vp — {m/2p)e~ m / 2 ' 5 , and the parameter 
K ~ (3 2 {p/m) 2 e m / 2 P eventually becomes large. Haldane's 
analyzis of the Bethe-ansatz equations has shown the 
crossover region at A(p/m) ln/3 _1 ~ 1 and the saturation 
of K. — * 1 at smaller densities. This statement, albeit 
reasonable, is hard to check in our formalism. 

The saturation of K can be a result of the interven- 
tion of some class of diagrams, increasingly important at 
small densities. As we have shown, the full Lagrangian 
( pq ) is free from infra-red divergencies. On the other 
hand, it requires the essential ultra-violet regularization. 
The situation is reversed when one works with the initial 



field variables ( |19[) . Hence the question of saturation of 
K. remains open within our formalism. 



VI. DISCUSSION AND CONCLUSIONS 



We have studied the sine-Gordon model in the pres- 
ence of a finite density of kinks. This quantum problem 
is related to the strongly interacting electronic system in 
one spatial dimension close to commensurability, where 
the metal-insulator transition occurs. The model was in- 
vestigated previously by the Bethe ansatz method, which 
gave a form of the spectrum and the Luttinger parame- 
ter, defining the decay of correlations functions at large 
distances. 

In dealing with complicated problems, it is always 
worth to have alternative methods at hand, suitable for 
cross-checking with existing results and capable for fur- 
ther theoretical predictions. In this paper we used the 
semiclassical approach and demonstrated that the spec- 
trum and the Luttinger parameter correctly reproduced 
the results of the Bethe ansatz studies. Among the new 
features, available due to the quasiclassics, are the ex- 
plicit form of various correlation functions and the un- 
ambigous construction of the charge-raising operator in 
the theory. 

Discussing the correlation functions, we determined 
their leading and subleading asymptotes, whose ampli- 
tudes and critical exponents strongly depend on the 
kink density. Addressing a more technical issue of the 
charge-raising operator, or the dual field in the theory, 
we showed the way it is constructed through the modi- 
fied collective coordinate method. This allowed us to re- 
establish the link with the fermionic counterpart of our 
bosonic problem and to calculate, in the leading order, 
the behavior of the superconducting- type correlations. 

One the consequences of our approach is the non-zero 
value of the quantum average, associated with the four- 
fcrmion umklapp operator in the Hubbard model away 
from half- filling. As a result, we found that close to com- 
mensurability point, the spectral density of the right- 
(left-) moving fermions is distributed over the neighbor- 
ing "4/ip" harmonics of the fermionic density. 
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APPENDIX: COLLECTIVE COORDINATE 
METHOD 

For brevity, we will write dtf = /, d x f = f in this 
section. Following p3[ we write our solution in the form 

oo 
4>(X, t) = (f>n(x + X) + ^ Cn(t)Vn{x + X) (Al) 



The classical solution <pn is proportional to /3 , and 
the fluctuations r\ are of order of (3°. In terms of the 
above mass of kinks condensate Mtot, the normalized 
zero-energy solution is 

r,o(x + X)=M; o ] /2 ci > ' Q (x + X), 

in accordance with (|37J). In order to pass to the correct 
form of the quantum Hamiltonian, one has to regard the 
variable X, which indicates the location of the soliton 
lattice, as a dynamical variable, replacing the zero mode 
coordinate cq. Instead of (Al) we write 



dC/diij — DjiUi. The classical Hamiltonian then 



cj>(x, t) = </>(,(£ + X{t)) + J2 Cn(t)Vn{x + X{t)), (A2) 

i 

and the time derivative 

oo oc 

j>(x, t) = ($, + ^2 CnV'n) X + ^ CnVn, (A3) 



1 



Introducing the new coordinates u n , n — 0, 1, . . . , oo 
according to u (t) = X(t), u n {t) — c n (t) for n > we 
represent the kinetic energy term in the Lagrangian in 

the form \ J Q dx <f> 2 — | X}i*7'=o UiD^jUj, with a symmet- 
ric matrix Da with elements 



becomes 
H 



2 t-^> 

i,3=0 



miD-^iM + VdUn}), 



(A5) 



where the potential energy term V is given by V({u n }) = 

/ dx[y 2 + um = M tot + \ EZi <f n + cm. 

The elements of the inverse matrix D 1 arc 



(D~ l ) 00 = l/D, 
(D-^on = -DoJD, n>0, 
(D~ 1 ) nm = S nm + D 0n D 0m /D, n, m > 0. 



(A6) 



and the determinant D of the matrix Dy is given by 
D = (y/Mtot + ao) 2 with a Q = J dxi]o(x)(J2^° u m r]' m ). 
Using above formulas, one obtains after some algebra 



<fi(x,t) 



Vo __ , v^/ , D 0r. 



7^0 



n=\ 



(Vr, 



Vo)n„ 



(A7) 



Requiring now [u n (t),ir m (t)] — iS nm and using 

7]), one verifies the exactness of the equation 
[<j)(x,t),<j)(y,t)] = i5(x — y), which result improves the 
previous formula (Kg). 

Among the new features, brought about by the de- 
scribed method is the explicit absence of the variable 
X(t) in the Hamiltonian ([Aq). As a result, its conjugate 
momentum ttq is independent of time and equals to the 
conserved total field momentum P = L dx 4><f>' , which 
equality is easily checked. 

If we consider the soliton lattice in its rest frame, 
P = 0, then the Hamiltonian is simplified and becomes 



n,m— 1 






)ir n ir m + V({«n}), (A8) 



At first glance, eqs. (A7), (A8) contain non-trivial ad- 
mixtures oc Don/y/T), which should interfere into the 
subsequent consideration. This is true, when one deals 
with a vacuum, consisting of a finite number of kinks Q. 
However, in the limit Q ~ L — > oo, one has r\ n oc L^ 1 ' 2 
and therefore D 0n = 0(1), whereas D = 0(Q(3~ 2 ). As 
a result, the terms D 0n /y/~D ~ (3Q- 1 / 2 in ([A?]), (|aJ) 
should be neglected in this limit and we arrive at sim- 
pler equations Ti = So + \ Y^=\ ( n n + u> n u n) + 0(fi) anc ^ 

4>(X, t) P=a = Y,n=l Vn(x + X)lT n . 



Dnn = dx 



1 



Dor, 



dx r) n C)^u m r]: 

l 
7l , n,m > 0. 



UnVn) 

mh n>0, 



(A4) 



The Hamiltonian is obtained by finding the canonical 
momenta ■Kj, conjugate to coordinates Uj, according to 
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